	do jj=1,8
	xt(1:3,jj)=meshPoints(1:3,meshCellPoints(jj,i))
	end do			
!	write(*,*) 'i=',i
	call meshCellOut_defn(xt,meshCellOut,nCellsOut,i,&
		areaOut(nCellsOut+1),xOut(1:3,nCellsOut+1),normalOut(1:3,nCellsOut+1))
